Zipf 's Law in the Liquid Gas Phase Transition of Nuclei 
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Zipf 's law in the field of linguistics is tested in the nuclear disassembly within the framework of 
isospin dependent lattice gas model. It is found that the average cluster charge (or mass) of rank n 
in the charge (or mass) list shows exactly inversely to its rank, i.e., there exists Zipf's law, at the 
phase transition temperature. This novel criterion shall be helpful to search the nuclear liquid gas 
phase transition experimentally and theoretically. In addition, the finite size scaling of the effective 
phase transition temperature at which the Zipf's law appears is studied for several systems with 
different mass and the critical exponents of v and /3 are tentatively extracted. 

PACS Number(s): 25.70.Pq Multifragment emission and correlations; 05.70.Jk Critical point phe- 
nomena; 24. 10. Pa Thermal and statistical models 



In intermediate energy heavy ion collisions (HIC) the 
hot nuclei with moderate temperature can be formed and 
finally de-excite by different decay modes, such as mul- 
tifragmentation. This kind of multifragment emission 
shows a rise and fall with beam energy (excitation en- 
ergy, or nuclear temperature, or impact parameter ...) 
HI , which probably relates to the critical point behavior 
in nuclear matter. The onset of multifragmentation, i.e., 
when the multiplicity Ni m f of intermediate mass frag- 
ment (IMF) attains to the maximum, probably indicates 
of the coexistence of liquid and gas phases. Relations to 
the IMF distribution, theory predicted that there exists 
a minimum of power law parameter, r TO j„, for the IMF 
cluster distribution when the liquid gas phase transition 
takes place Q . 

On the other hand, the study on nuclear caloric curve 
is an attractive subject due to the fact that a tempera- 
ture plateau exists when the first order liquid gas phase 
transition occurs. Experimentally He-Li isotopic temper- 
ature from Albergo's thermometer || for projectile-like 
Au spectators seems to have a temperature plateau in the 
excitation energy range of 3 to 10 MeV/u by ALADIN 
Collaborations M , which was claimed that the first order 
nuclear liquid gas phase transition was observed exper- 
imentally. Nuclear caloric curves were also surveyed by 
several groups, such as INDRA [§, MSU @, TAMU § 
, EOS § and GSI § etc. 

In addition, the possibility to extract critical exponents 
and critical behavior from finite-size systems has been 
attempted experimentally in [ 10|] - 

However, even though the above trials to search for the 
liquid gas phase transition are multi- variant, these signa- 
tures are still in controversial debating JlT| p^|. In this 
context, it will be helpful and meaningful to present novel 
signature to characterize the nuclear liquid gas phase 
transition to guide the experimental analysis and the- 
oretical predictions. 

In this paper we will introduce, for the first time, Zipf's 
law into the study of nuclear liquid gas phase transi- 
tion. Zipf's law pa has been known as a statistical phe- 



nomenon concerning the relation between English words 
and their frequency in literature in the field of linguis- 
tics. In this Zipf's analysis, one calculates the histogram 
that gives the total number of occurrences of each word 
in a text. If all the words in the text are arranged in 
rank n order, from the most frequent to the least fre- 
quent, then such a histogram is found to be linear on 
double logarithmic paper, with a slope -A, with A f» 1 for 
all languages. Attempts to understand the origin of the 
Zipf's law are connected to the hierarchical structure of 
language |T(| . This relation was found not only in linguis- 
tics but also in other fields of sciences. For examples, the 
law appeared in distributions of populations in cities, dis- 
tributions of areas of lakes, DNA base pair sequences |17|| 
and cluster-size distribution in percolation process p8| . 
In this paper, Zipf's law will be tested for the charge (or 
mass) distribution of nuclear clusters and evidenced to 
be a factor to characterize the phase transition. 

Isospin dependent lattice gas model (I-LGM) was used 
to investigate the Zipf's law in nuclear cluster distribu- 
tion. The lattice gas model was developed to describe 
the liquid- gas phase transition for atomic system by Lee 
and Yang |1!| first. The same model has already been 
applied to nuclear physics for isospin symmetrical sys- 
tems in the grancanonical ensemble pp|] with a sampling 
of the canonical ensemble |ll],^l]-^6| , and also for isospin 
asymmetrical nuclear matter in the mean field approxi- 
mation p7{ . Here we will make a brief description for the 
I-LGM model. 

In the isospin dependent lattice gas model, A (= N+Z) 
nucleons with an occupation number s which is defined 
s = 1 (-1) for a proton (neutron) or s — for a vacancy, 
are placed on the L sites of lattice. Nucleons in the near- 
est neighboring sites interact with an energy e SiSj . The 
hamiltonian is written as 
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The interaction constant e SiSj is chosen to be isospin de- 
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pendent and be fixed to reproduce the binding energy of 
the nuclei 25 1: 



< pp 
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which indicates the repulsion between the nearest neigh- 
boring like-nucleons and attraction between the nearest 
neighboring unlike- nucleons. This kind of isospin depen- 
dent interaction incorporates, to a certain extent, Pauli 
exclusion principle and effectively avoids producing un- 
reasonable clusters, such as di-proton and di-neutron etc. 
Three-dimension cubic lattice with L 3 sites is used which 
results in pj = j^po of an assumed freeze-out density of 
disassembling system, in which po is the normal nuclear 
density. The disassembly of the system is to be calcu- 
lated at pf, beyond which nucleons are too far apart 
to interact. Nucleons are put into lattice by Monte 
Carlo Metropolis sampling. Once the nucleons have been 
placed we also ascribe to each of them a momentum by 
Monte Carlo samplings of Maxwcll-Boltzmann distribu- 
tion. 

Once this is done the I-LGM immediately gives the 
cluster distribution using the rule that two nucleons are 
part of the same cluster if 



P 2 /2p 



< 0. 
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This prescription is evidenced to be similar to the 
Coniglio-Klein's prescription (2^| in condensed matter 
physics and shown to be valid in lattice gas type cal- 
culations iinim]. 

We first chose the medium size nuclei 129 Xe to analyze 
the phase transition point behavior in nuclear disassem- 
bly and its Zipf's law from the cluster distribution. The 
freeze-out density pf is chosen to be 0.38 po due to the 
data were best fitted by a pf between 0.3po and 0.4po 
in the previous LGM calculations [ p3p9| ], which corre- 
sponds to the cubic lattice with size L = 7 for 129 Xe. 
f 000 events are simulated for each calculation. 

In order to check the phase transition behavior in the 
I-LGM, we will firstly show the results of some physi- 
cal observables, namely the effective power-law parame- 
ter, r, the second moment of the cluster distribution, 52 
J30J, and the multiplicity of intermediate mass fragments, 
Ni m f for the disassembly of 129 Xe in figure 1. These ob- 
servables have been evidenced useful in previous works 
to judge the liquid gas phase transition, as shown in Ref. 
p5| , pl|j32| . The valley of r, the peaks of N irn f and 52 
happens around T ~ 5.5 MeV which is the signature of 
onset of phase transition. However, the aim of this paper 
is searching a novel signature of liquid gas phase tran- 
sition besides the above observables. The above phase 
transition temperature will be only used as a reference of 
the novel signature, as stated below. 

Now we present the results for testing Zipf's law in the 
charge distribution of clusters. The law states that the 
relation between the sizes and their ranks is described by 



Z n = c/n (n=l, 2, 3, ...), where c is a constant and Z n 
is the average charge (or mass) of rank n in a charge (or 
mass) list when we arrange the clusters in the order of 
decreasing size. For instance the charge Zi of the second 
largest cluster with rank n — 2 is one-half of the charge 
Z\ of the largest cluster, the charge Zz of the third largest 
cluster with rank n = 3 is one-third of the charge Z\ of the 
largest cluster, and so on. In the simulations of this work, 
we averaged the charges for each rank in charge lists of 
the events: we averaged the charges for the largest clus- 
ters in each event, averaged them for the second largest 
clusters, averaged them for the third largest clusters, and 
so on. From the charges averaged, we examined the re- 
lation between the charges Z n and their ranks n. Figure 
2 shows such relations of Z n and n for Xe in different 
temperature. The histogram is the simulated results and 
the straight lines represent the fit with Z n oc n~ A in the 
range of 1 < n < 10, where A is the slope parameter. A is 
5.77 at T = 3 MeV. Then we increased the temperature 
and examined the same relation and obtained A = 3.65 
and 1.53 at T = 4 and 5 MeV, respectively. Up to T = 
5.5 MeV, A = 1.00, i.e., at this temperature the relation 
is satisfied to the Zipf's law: Z n oc n . When tempera- 
ture continues to increases, A continues to decreases, for 
instance, A = 0.80 at T = 6 MeV and A = 0.56 at T = 
7. This temperature having the Zipf's law, denoted as 
Ta, is consistent with the phase transition temperature 
obtained in Fig. 1, illustrating that the Zipf's law is also 
a good judgement to phase transition. From the statisti- 
cal point of view, the Zipf's law is related to the critical 
phenomenon PJ33]] . Figure 3a summarizes the parameter 
A as a function of temperature. Clearly the Zipf's law (A 
= 1) reveals at phase transition point. 

In order to further illustrate that the Zipf's law exists 
most probably in phase transition point, we directly re- 
produce the histograms with Zipf's law: Z n = c/n. In 
this case, c is sole parameter, but what we are interest- 
ing in is its truth of the hypothesis of Zipf's law: the \ 2 
test. Figure 3b demonstrates the \ 2 l n df for the Z n - n 
relations at different T. As expected, there is the min- 
imum x 2 / n df around the phase transition temperature, 
which further support that Zipf's law of the fragment 
distribution reveals when the liquid gas phase transition 
occurs. 

In the above examinations, the system size is finite, so 
the effective phase transition temperature, Ta, deduced 
from the Zipf's law shall be only valid for a specific fi- 
nite size nuclei. As well known, the finite size scaling is 
an important factor in studying the phase transition of 
finite matter. In this context, we will investigate such 
an effect for Ta with the help of the Zipf-type analysis 
for nuclear clusters and predict the value of critical tem- 
perature, Tc, at which the Zipf's law will appear for the 
infinite nuclear matter with the same N/Z and freeze-out 
density as 129 Xe in the framework of isospin dependent 
lattice gas model. When we examine the values Tc of 
thresholds of phase transition, we usually extrapolate the 
effective phase transition temperature Ta as L — > oo [pil . 
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This is based on the extrapolation rule 

\T C -T A \ ocL^, (4) 

where v is the critical exponent for correlation length. 
For three dimension (3D) percolation class, v = 0.9 ; for 
the 3D Ising class or liquid gas class, v = 0.6 p5| . Below 
we will examine if the finite size scaling law is valid for 
the phase transition temperature at which the Zipf 's law 
occurs in the nuclear system and tentatively extract the 
critical exponents v and (3. 

We simulated the nuclear disassembly for systems with 
the same N/Z as 129 Xe but A = 80, 274, 500, 830 and 
1270 within the cubic sites of 6 3 , 9 3 , ll 3 , 13 3 and 15 3 , re- 
spectively, which will result in the same freeze-out density 
as 129 Xe for comparison. The effective phase transition 
temperatures, Ta, at which Zipf's law (A = 1) appears 
are extracted for each system. The figure 4a summarizes 
such Ta as a function of L, from which we can get the 
two fit parameters of Tc and v with the Eq.(4). The 
best fit parameter is shown in the figure, i.e., Tc = 5.846 
± 0.013 MeV and v = 0.348 ± 0.023. It tells us that 
the scaling law (Eq.(4)) is rather valid even for the ef- 
fective phase transition temperature deduced from the 
Zipf's law. Viewing from the fit values of v, it seems to 
approach to the value of 3D Ising or liquid gas class than 
the 3D percolation class, which probably indicate that 
the nuclear disassembly belongs to 3D liquid gas univer- 
sality class, while the corresponding extrapolated critical 
temperature of Tc — 5.846 MeV for the infinite nuclear 
matter with N/Z = 1.39 and 0.38po as A — > oo. In ad- 
dition, the finite size scaling can be further investigated 
by the source size dependent largest cluster at the phase 
transition point where the Zipf's law appears and hence 
to extract the critical exponent (3 by the scaling law: 

A max /L 3 oc L-V", (5) 

where A max is the size of the largest cluster at the phase 
transition point where A = 1 and L is the lattice size 
of studied system. For three dimension (3D) percola- 
tion class, (3/v = 0.44; for the Ising class or liquid gas 
class, (3/v = 0.51 pj|. Figure 4b shows the scaling law 
of Eq.(5) for A = 80, 129, 274, 500, 830 and 1270. This 
method leads to (3/v = 1.105 ± 0.039 which is close to 
the classical value. However, we would like to emphasize 
that one should not be nutty to the value itself of the 
exponents because they are very sensitive to the reduced 
Ta and A max . In addition, the size of lattice should be 
as large as possible if one really wants to extract the pre- 
cise exponents. What we want to say here is that the 
finite size scaling is valid for the Zipf-type analysis and 
it is also possible, in principal, to study the critical ex- 
ponents from there. 

In conclusion, Zipf-type analysis in the field of linguis- 
tics is applied to the cluster charge (mass) distribution 
which has been generated with an isospin dependent lat- 
tice gas model. The cluster distributions at the phase 
transition point show exactly inversely to its rank, i.e., 



there exists Zipf's law. This criterion shall be useful in 
searching the nuclear liquid gas phase transition experi- 
mentally and theoretically in combination with the other 
phase transition observables, such as t, S2 and N m f etc. 
In addition, the finite size effect of the effective phase 
transition temperature at which the Zipf's law appears 
is also addressed and the critical exponents of v and (3 
are tentatively extracted. It seems that the nuclear dis- 
assembly approaches to the 3D liquid gas universality 
class in the framework of the isospin dependent lattice 
gas model. 
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Fig. 4: The effective phase transition temperature (a) 
and the normalized largest cluster size A max /L 3 (b) is 
plotted as a function of size of lattice. From (a) the crit- 
ical temperature Tc for infinite system and the critical 
exponent v can be extracted; From (b) the critical expo- 
nent (5/v can be obtained. 



Fig.l: The effective power-law parameter, r, the sec- 
ond moment of the cluster distribution, S2 , and the mul- 
tiplicity of intermediate mass fragments, Ni m f as a func- 
tion of temperature for the disassembly of 129 Xe. The 
arrow represents the position of phase transition temper- 
ature. 

Fig. 2: The average charge Z n with rank n as a function 
of n for 129 Xe. The histograms are the calculation results 
and the straight lines are their fits with Z n oc n 

Fig. 3: The slope parameter A of Z n to n (a) and the 
X 2 test for Zipf's law (b) as a function of temperature 
for 129 Xe. The arrow represents the position of phase 
transition temperature. 
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